AD-A206  670 


i 


OFFICE  OF  NAVAL  RESEARCH 
Contract  N00014-86-K-0043 
TECHNICAL  REPORT  No.  96 

Light  Scattering  from  an  Atom  Near  the  Surface  of  a  Superlattice 
by 

Xiao-shen  Li,  D.  L.  Lin  and  Thomas  F.  George 

Prepared  for  Publication 
in 

Atomic  and  Molecular  Physics 

Edited  by  M.  S.  Z.  Chagtai 

Aligarh  Muslim  University,  Aligarh,  India 


Departments  of  Chemistry  and  Physics 
State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


April  1989 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


dtic 

electe 

APR  1  1 198S 


89  4  10  016 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 

UBUFFALO/DC/ 89/TR-96 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No.  0704-0188 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Depts.  Chemistry  &  Physics 
State  University  of  New  York 


6c  ADDRESS  (City,  State,  and  ZIP  Code) 

Fronczak  Hall,  Amherst  Campus 
Buffalo,  New  York  14260 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


6b.  OFFICE  SYMBOL 
(If  applicable) 


8b.  OFFICE  SYMBOL 
(If  applicable) 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 

unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 
Chemistry  Program 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Contract  N00014-86-K-0043 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


Office  of  Naval  Research 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 
Chemistry  Program 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


1 1 .  TITLE  (Include  Security  Classification) 

Light  Scattering  from  an  Atom  Near  the  Surface  of  a  Superlattice 


PROJECT 

TASK 

NO. 

NO 

WORK  UNIT 


12.  PERSONAL  AUTHOR(S) 


13a.  TYPE  OF  REPORT 


Xiao-shen  Li,  D.  L.  Lin  and  Thomas  F.  George 


13b.  TIME  COVERED  114.  DATE  OF  REPORT  (Year,  Month,  Day)  |15.  PAGE  COUNT 

FROM _ TO _  I  I 


16.  SUPPLEMENTARY  notation  Prepared  for  publication  in  Atomic  and  Molecular  Physics,  Edited  by 

M.  S.  Z.  Chagtai,  Aligarh  Muslim  Univeristy,  Aligarh,  India 


COSATI  CODES 


GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

ATOM,,  ^  LIGHT  SCATTERING/ 

NEAR4SUPERLATTICE  RESONANCE  FLUORESCENCE 

METAL-INSULATOR  _ -  '  SEMICLASSICAL  THEORY. 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

'""A  semiclassical  theory  is  developed  to  treat  the  interaction  of  radiation 
with  an  atom  located  near  the  surface  of  an  m-component  superlattice .  For  the 
special  case  of  m  =  2,  the  radiation  penetration  depth  into  the  superlattice 
and  the  resonance  fluorescence  spectrum  are  calculated  numerically  for 
superlattices  of  metal-insulator  construction.  The  results  show  sensitive 
dependence  on  the  plasma  frequency  in  the  metal  layer  as  we] 1  as  on  the 
dielectric  function  and  width  of  the  insulator  layer.  ^  r 


20.  DISTRIBUTION  /  AVAILA8ILITY  OF  ABSTRACT 
SI  UNCLASSIFIED/UNLIMITED  0  SAME  AS  RPT 


22«.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  David  L.  Nelsoi 


DO  Form  1473,  JUN  86 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  users  Unclassified 


22b.  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 
( 2021  6Q6-441D 


Previous  editions  are  obsolete.  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


I.  Introduction 


When  an  atom  or  molecule  near  or  adsorbed  on  a  solid  surface  interacts 
with  a  radiation  field,  the  presence  of  the  surface  and  substrate  greatly 

influences  the  optical  phenomena.  Examples  include  surf  ace  -  enhanced  Raman 

12  .  .  3-6 

scattering,  ’  coherence  and  energy  transfer  in  spontaneous  emission, 
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surface- induced  resonance  fluorescence,  and  so  on.  The  optical  properties 

of  adspecies  on  a  solid  surface  may  provide  a  sensitive  probe  of  electronic 
and  other  structure  of  the  substrate,  and  hence  has  prompted  a  careful  re¬ 
examination  of  the  optics  of  surfaces.  Once  the  mechanism  of  these  optical 
processes  is  fully  understood  and  brought  under  experimental  control,  it 
should  become  a  powerful  tool  for  the  analysis  of  surface  processes. 

On  the  other  hand,  when  a  strong  driving  coherent  field  is  nearly  on 
resonance  with  the  adatomic  transition  frequency,  the  field  places  the  driven 
atom  in  an  environment  in  which  the  probability  of  stimulated  emission  may 
exceed  that  of  spontaneous  emission.  When  this  is  the  case,  ac  Stark 
splitting  of  the  atomic  levels,  Rabi  oscillations  of  the  level  occupation 

probabilities  and  nutational  oscillations  of  the  stimulated  field  intensity 

7-14 


are  enhanced.  Consequently,  resonance  fluorescence 

15,16 


and  other  nonlinear 

optical  phenomena  become  possible."”'’"''  Because  multiphoton  transitions  are 
just  as  important  as  the  s ingle  -  photon  transition  in  such  light- driving 


processes,  the  ordinary  perturbation  method  is  no  longer  reliable. 

In  the  surface -free  case  such  as  atoms  in  a  gas , ^  such  processes  can  be 
effectively  dealth  with  by  the  well-known  optical  Bloch  equation  (OBE) .  With 
the  presence  of  the  surface,  a  set  of  surface-dressed  optical  Bloch  equations 
(SBE)  has  been  derived  ^  ^  to  treat  such  problems  as  the  effects  of  surface- 
reflected  photons,  the  resonance  interaction  of  an  adatom  with  surface 
plasmons,  collisional  dephasing  of  the  adatom  due  to  the  gas  atoms  in  the 


medium,  and  random  phase  fluctuation  of  the  intense  laser  field.  Resonance 
fluorescence  has  been  investigated  for  an  adatom  near  a  flat  metallic 


7  S 

surface  J  and  near  a  rough  metallic  surface  represented  by  a  hemispheroid 

rA  1 A 
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protrusion  on  a  flat  surface 


More  recently,  a  different  approach 


involving  reservoir  theory  and  Dekkars ' s  quantization  procedure  for  a 
dissipative  system'1^’^  has  been  adopted  to  derive  another  set  of  the  SBE.  " 
The  resonance  fluorescence  spectrum  of  an  adatom  has  been  analyzed  for  the 
following  cases:  an  adatom  is  considered  near  the  surface  of  a  metal  sphere 

n 

with  the  size  of  the  sphere  taken  into  account,  and  the  problem  is  considered 

again  near  the  surface  of  a  semi-infinite  conductor  and  dielectric  superlattice 

in  which  the  frequency  of  the  transverse  optical  phonon  matches  that  of  the 
.  .  11 

adatomic  transition. 

We  investigate,  in  this  paper,  the  resonance  fluorescence  spectrum  of  an 

adatom  near  the  surface  of  a  superlattice.  The  adatom  is  taken  as  an  emitting 

dipole.  The  emitted  field  is  reflected  back  from  the  surface  and  interfaces 

and  is  coupled  to  the  radiating  dipole  whose  dynamical  behavior  is  therefore 

totally  changed.  We  first  develop  the  theory  for  a  general  superlattice 

composed  of  m  different  layers  per  period.  This  theory  may  be  applied  to 

20  21 

simulate  quasi-periodic  systems.  ’  Numerical  results  are  given  for  the 
special  case  of  m  -  2. 


II.  Theory 

For  a  two- level  atom  driven  by  a  monochromatic  laser  field  with 


E(t)  -  ~  (Ee 


+  E*e 


(1) 


4 


the  SBE  are  given  by 
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[<sV|  fi(A+nS)  -  7  10 


f  <sz> 

dt 

<S> 


i0*/2 


■2y  - i0/2 


f<S+>l 


<s2> 


-iO*  - i(A+0  )  -  7 1  L <S  >1 


The  notation  is  as  follows.  The  adatom  with  a  transition  frequency  w  is 
located  at  a  distance  dQ  away  from  the  surface  of  the  super lattice .  The 
matrix  element  of  the  electric  dipole  moment  operator  is  denoted  by  |p|,  and  E 
and  are  the  amplitude  and  frequency  of  the  external  laser  field, 
respectively.  The  detuning  is  A  -  u  -  u>  ,  and  the  Rabi  freuquency  is  Q  - 
| p | E .  The  transition  probability  amplitude  is  proportional  to  the  projection 
operators  defined  by 


S  -  l+X- 


$  (|+x+j 


-X-  ) 


The  total  decay  rate  of  the  adatom  can  be  written  as 


o  s 
7-7  +7 


o  22 

where  7  is  the  decay  rate  in  the  absence  of  the  substrate, 


o  ^  y —  I  1 2  3  .  3 
7  “  3  'P'  w  /C 
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and 


yS  =  |p| 2  Im  f(dQ) 


(4b) 


is  the  decay  rate  induced  by  the  surface.  The  frequency  shift  of  the 
spontaneous  radiation  due  to  the  surface  is 


nS  -  I P I  2  Re  f(do)  , 


and  the  function  f(d  )  is  determined  bv  ^  ^ 

o 


(5) 


Er  =  | P | f (dQ) S "  =  pf(dQ)  ,  (6) 

where  E  is  the  component  of  the  reflected  field  E  in  the  direction  of  p. 

K  R 

Equation  (2)  agrees  with  the  result  from  linear  response  theory^  when 

the  adatom  is  taken  as  a  harmonic  oscillator.  With  the  aid  of  the  regression 

23  24- 

theorem  for  correlation  functions,  5  we  find  from  the  SBE  the  well-known 

H-IH 

results  of  the  incoherent  resonance  fluorescence  spectrum, 

^inc(l/)  "  2  lnl^(I)2  +  2  |n|2  +  ^Tf2>/<2  l°l 2  +  |z|2>  (x2+y2)  ,  (7) 

12  2  2 

where  D  -  u  -  ,  z  -  7  +  i(A+nS),  x  -  2  -7  ( ~  |  CJ  j  +  |z|  -  2D  ),  and 

y  -  D(|0|2  +  | z | 2  +  472  -  D2). 

Consider  now  a  semi - inf inite  superlattice  consisting  of  m  different 
layers  per  period  as  shown  in  Fig.  1.  The  surface  is  taken  to  be  the  xy- 
plane,  and  an  adatom  with  dipole  moment  p  is  located  at  r  -  r  ■  (0,0, -d^). 
The  dielectric  constants  and  thicknesses  of  these  layers  are  and  d^ , 
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respectively,  where  i  -  1,  2,  .  .  . ,  m.  The  reflected  electric  field  at  rQ  from 
a  substrate  of  a  two-compo;i’nt  superlattice  has  been  worked  out  in  Ref.  8.  We 
now  generalize  the  theory  to  the  case  of  an  m-component  superlattice  in  which 

■4 

the  thickness  of  each  period  is  L.  The  set  of  equations  for  the  fields  E  and 

H  in  different  materials  can  be  obtained  directly  from  Maxwell's  equations. 

i  G, 11 .  13  14 

The  results  are  ’  -  5 


2-»  2-*- 
V  E  +  e  k  E 
o  o 


2->  1 

-  Un[k  P  +  “  7(V*P) 

O  £ 

O 


H  =  TT~  V  x  E  , 
ik 


z  <  0 


(8a) 


2-*  1-* 

VZE.  +  £  k  E,  =  0 
1  1  o  1 


V  •  Ex  -  0 


H1  ’  ik“  V  x  E1  ’ 
o 


nL  <  z  <  nL  +  d. 


(8b) 


2-»  2-* 
v  E_  +  £.k  E.  =  0  , 

Z  Z  O  Z 


V  •  E2  -  0 


H2  -  Ik^  7  X  E2 


nL  +  d^  <  z  <  nL  +  d^  +  d^ 


(8c) 


2->  2-* 
v  E  +  £  k  E  -  0  , 

m  mom 


V  •  E  -  0 
m 


H  -  -f-  7  x  E 
m  ik  m 

o 


(n+l)L  -  dm  <  z  <  (n+l)L  , 


(8d) 


where  kQ  -  w/c  is  the  wave  number  in  vacuum,  n  labels  the  period  in  the 

4 

superlattice,  and  P  is  defined  by 


(9) 


P(r,w)  -  p(w)  £(r-rQ) 


The  solution  for  the  electric  field  in  Eq .  (8a)  can  be  written  as 


j  w)  -  J  du  I  dv  £(u,v,cd)  e^k  r  +  E^Cr.w)  ,  z  <  0 


-  0  , 


k'  -  (k«,-w) 


2  ,  ,2  ,2 
w  -  €  k  -  ki| 
o  o 


ImW  >  0 


where  we  have  defined  k  -  (kii ,w)  -  (u,  v,  w) .  The  field  due  to  the  dipole  is 


E  (r,w) 
P 


duI 


2-*  _  - 1  -» 
k  p  -  £  V(p*7) 

o  o _ 


x  exp[iu(x-XQ)  +  iv(y-yQ)  +  iw|z-ZQ|] 


Because  of  the  periodic  structure  of  semi -infinite  superlattice,  the  electric 
field  may  have  a  Bloch-wave-like  form  with  an  envelope  function  decaying 
exponentially  with  increasing  z.  Thus  we  find  solutions  for  the  electric 
field  of  Eqs .  (8b) -(8d)^  ’  3S 


iw. ( z - nL) 


Ef(r  ,w)  -  |  du  J  dv  (u,v,u)e  1 


-*■  -iw  (z-nLK  i£i.*r  -  /9nL 

+  E^Oi.v.w)  e  1  }e  « 


-♦  Q  (+)  r*  C  (  \ 

ki  ‘  £l  °  4  *  ^  s 


(12a) 
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kl  -  (kii  ,w1)  ,  k£  -  (k|. ,  -wi) 


2  2  2 

W1  “  elko  ‘ 


nL  <  z  <  nL  +  d, 


E2(r,w)  -  J  du  J  dv  j^2+^(u,v,w) 


iw2(z-nL-d1) 


+  \u,v,u)e  2 


iw0(z-nL-d.  )n  ik|.*r  -  0nL 


E2.£^-E2.£^-o 


k2  “  <k||  ’  k2  "  (k|!  ’  "'V 


2  2  2 
"2-‘2ko-k  • 


nL  +  d^  <  z  <  nL  +  d^  +  d2 


(12b) 


Em(r,w)  -  J  du  J  dv  {^+)(u,v,U) 


iw  [z  -  (n+l)L  +  d  ] 


7?,  .  -lwm[ z  -  (n+l)L  +  d  k  ik(|*r  -  /3nL 

+  £<-><«, v,„)  e  “  ■  }  e  I 

m  J 


k  .  £  (+)  -  k  •  £  (_)  -  o 

mm  mm 


k  -  (kii  ,W  )  k'  -  (kii  ,  -w  ) 

m  m  m  m 


2  i2  »  2 

w  -  e  k  -  Kii 
m  mo  i! 


(n+l)L  -  d  <  z  <  (n+l)L 
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(12c) 


where  we  have  introduced  a  parameter  p  to  measure  the  rate  of  absorption. 
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For  simplicity,  we  assume  that  the  dipole  is  orientated  perpendicular  to 
the  interfaces  of  the  super lattice .  We  can  then  follow  the  standard  procedure 
to  find  the  component  of  the  reflected  field  along  the  dipole  direction  at  r  . 
Equations  (12)  together  with  the  usual  Maxwell  boundary  conditions  at  the 
surface  and  interfaces  then  yield 


E 


R 


r  ko  p  0(do) 


(13) 


0(do) 


2  iU  d 
o  o 


e 


V*2 


(14a) 


-  *  2m'1  (1  -  V1Q)  e^L 


1  {■ 


N 


+  )  1(1*  (-D  ra  V  j 

L  mo 

N  -0 
m 


y  {expu}  (-i>  1  ujdj 

N2’*"’Nm-l'°  j"1 


VN2+i 

x  n  [l  +  (-D  v 


m- 1 

n 

2=1 


2,2+1 


■}} 


(14b) 


*2  ’  2”'l<1  *  V10> 


1  N  1  m  N 

}  {[i  *  (.i,  \0i  }  f i  \  /  -  1  \  J 


N  =0 
m 


exp[ i  )  (-1)  J  U  d ] 


VA-l"0  J"1 


ra'1  ri  Nf+Nf+i  n 

W  (  >  vi,i+i]}}  - 


(14c) 


where  N]_  =  =  0,  U,  -  Je £- k  ,  d,  -  kd,,  and 
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If  the  superlattice  is  composed  of  highly-dissipative  materials,  the  incident 
electric  field  cannot  go  very  deep  into  the  layers.  When  the  condition  L  «  X 
is  satisfied,  which  is  usually  the  case,  we  can  define  the  penetration  depth 
for  the  field  as 

A 

dp  -  aL  -  (Re$)  ,  (17) 
where  a  is  a  positive  number. 

It  is  observed  from  the  above  equations  that  both  the  reflected  field 
and  the  penetration  depth  depend  strongly  upon  the  number  of  layers  per  period 
as  well  as  the  thickness  and  dielectric  properties  of  the  individual  layer. 
Meanwhile,  the  surface  plasmon  TM  modes  are  determined  by  the  zeroes  of  the 
denominator,  namely,  -  0.  The  properties  of  the  surface  plasmon  as  well  as 
the  interference  between  the  electric  fields  reflected  and  refracted  at  each 
interface  play  very  important  roles  in  the  determination  of  E^,  both  its 
amplitude  and  its  phase. 
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To  find  the  spontaneous  decay  rate  and  frequency  shift,  we  first  note 

that  a  comparison  of  (6)  and  (13)  yields  f(d  )  -  kJ  <p(d  )•  Substituting 

0  6  0  0 
O 

into  (4)  and  (5),  we  obtain 

7  -  1  +  J  £*V2  Re  *(d  )  (18a) 

2  o  o 

nS  -  -  J  e'3/2  Im  4>(d  )  (18b) 

2  o  o 

expressed  in  the  unit  of  7°. 

III.  Numerical  results 

Th  •.  theory  developed  above  is  valid  for  any  number  of  layers  in  one 

period  of  the  superlattice.  In  what  follows,  we  shall  apply  it  to  a  typical 

case  of  two-component  superlattices .  The  first  layer  in  each  period  is  a 

2  2 

metal  with  dielectric  function  (v)  -  1  -  w  /(w  +  iu>7  )  ,  while  the  other  is 

1  p'  p 

an  insulator  with  a  wide  energy  gap.  Setting  m  =»  2  in  Eqs .  (14b),  (14c)  and 
( 16) ,  we  have 


or  iU  d 

+!  -  -2(l-v10)  e  +  [(1-V20)(l+V12)e 


A  A 


-iu2d2  iUxd 

+  (1+V20)(1-V12)  e  He  1 


m  iU2d? 

-  2(l+V10)ePL  -  [(1+V20)(l+V12)e 


(19a) 


(19b) 


n 


+  (1-V20)(l-V12)e  ']  e 

-l  A  A 

-  j  (l+V12)(l+V21)cos(U1d1  +  U2d2) 

A  A 

+  (l-V12)(l-V21)cos(U1d1-U2d2) 


(20) 


We  look  at  the  penetration  depth  for  normal  incidence  of  light  with 

A 

frequency  w.  From  (15),  (17)  and  (20),  we  calculate  d  as  a  function  of  the 

P 

adatomic  transition  frequency.  The  results  for  different  layer  thicknesses 
are  shown  in  Fig.  2,  in  which  we  have  defined  x  —  .  It  is  observed  that 

for  fixed  u,  the  penetration  depth  has  a  sensitive  dependence  on  the  charge 
carrier  density  in  the  metal  layers.  It  decreases  with  increasing  carrier 
density,  which  is  easy  to  understand  because  a  large  carrier  density  results 
in  more  collisions  and  therefore  more  energy  loss  of  the  incident  radiation. 
Furthermore,  we  also  see  that  the  larger  the  fraction  that  the  nonabsorptive 
insulator  in  each  period  occupies,  the  deeper  the  field  penetrates  into  the 
superlattice . 

The  incoherent  part  of  the  resonance  fluorescence  spectrum  may  be 
investigated  numerically.  We  find  again  that  the  spectrum  is  extremely 
sensitive  to  the  plasma  frequency  or  the  density  of  charge  carriers  in  the 
metallic  layers.  When  w  is  fixed,  different  carrier  densities  may  result  in 
spectra  with  different  number  and  shape  (height  and  width)  of  peaks  as  well 
as  the  position  of  sidebands.  The  situation  is  illustrated  in  Fig.  3  in  which 

the  spectrum  is  shown  as  a  function  of  D  [defined  just  beneath  Eq.  (7)]  for 
various  values  of  x.  That  the  sidebands  do  not  appear  in  the  case  of  the 
one-peak  spectrum  may  be  understood  as  follows. 
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When  x  is  such  Chat  4>0  is  near  its  minimum,  the  resonance  condition  for 
the  surface  plasmon  to  exist  is  nearly  satisfied.  The  creation  of  the  surface 
plasmon  then  enhances  the  absorption  and  results  in  a  single-peak  spectrum. 
This  is  because  non- radiative  energy  transfer  to  the  surface  plasmon  becomes 
very  strong  near  resonance  and  hence  greatly  increases  the  absorption  of  the 
incident  energy  by  the  surface  via  the  adatom.  Consequently,  the  scattered 
field  intensity  is  too  weak  to  produce  the  sidebands  of  the  spectrum.  On  the 
other  hand,  when  x  is  such  that  ^  is  far  from  its  minimum,  this  absorption 
process  is  insignificant,  and  the  spectrum  has  three  peaks  which  are  all  high 
and  narrow,  provided  that  there  is  no  significant  absorption  resulting  from 
the  multibeam  interference. 
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2.0,  |n|  -  45,  7p  -  0.01  M  £2  -  3.0,  -  0.1,  d3  -  0.5  and  (a)  x  -  0.1, 

(b)  x  -  0.2,  (c)  x  -  0.3,  (d)  x  -  0.4. 
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